A combinatorial solution for the current fluctuations in the exclusion process 
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We conjecture an exact expression for the large deviation function of the stationary state current 
in the partially asymmetric exclusion process with periodic boundary conditions. This expression is 
checked for small systems using functional Bethe Ansatz. It generalizes a previous result by Derrida 
and Lebowitz for the totally asymmetric exclusion process, and gives the known values for the three 
first cumulants of the current in the partially asymmetric model. Our result is written in terms of 
tree structures and provides a new example of a link between integrable models and combinatorics. 
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I. INTRODUCTION 

The asymmetric simple exclusion process (ASEP) is one of the simplest interacting particles systems featuring an 
out of equilibrium stationary state. It has been studied much in the past 0,0,0], in particular because it belongs to the 
class of exactly solvable models. A quantity of interest is the macroscopic stationary state current and its fluctuations, 
since the presence of this current is the signature that the system is out of equilibrium. Various boundary conditions 
have been used in the study of the one-dimensional ASEP: open boundaries connecting the system to reservoirs of 
particles 0, H, Q , infinite line Z 0, H, Q , and periodic boundary conditions [13, flll.fT^ , which is the case studied here. 

Using the Bethe Ansatz, all the cumulants of the current were calculated [ij] in the special case of the totally 
asymmetric simple exclusion process (TASEP), for which the particles hop in only one direction. For an arbitrary 
asymmetry between the hopping rates, finite size expressions for the three first cumulants were derived [Tst. and 
all the cumulants were obtained in the large system size limit with non vanishing asymmetry. 

In the present work, we conjecture an expression (|27|) generalizing these results: it provides an exact expression 
for all the cumulants of the stationary state current in the ASEP with partial asymmetry for finite systems. This 
expression gives the TASEP result [l^ and the limit obtained in [l^ . It also allows the study of a vanishing asymmetry 
probing the transition between the equilibrium system with symmetric rates and the totally asymmetric system for 
which detailed balance is maximally broken. The exact formulas obtained previously for the three first cumulants 
in the partially asymmetric model are also recovered. Our conjecture is checked for small systems using functional 
Bethe Ansatz. 

We consider the asymmetric simple exclusion process on a ring of size L with n particles hopping locally both one 
site forward (with rate p) and backward (with rate q = xp). By the exclusion rule, the particles are only allowed 
to hop if the destination site is empty. The integrated current Yt is defined as the total distance covered by all the 
particles between time and time t. When t becomes large, the fluctuations of Yt are given by [3, 

(e-/^') ^e^(7)* . (1) 
The formal series i?(7) is the exponential generating function of the cumulants of the current: 

E{l) = Jl+^l^ + ^l^ + ..., (2) 

where J is the mean value of the current, D the diffusion constant and Ey, the third cumulant of the current. The 
generating function E['^), which is related to the large deviation function of the current by a Legendre transform, is 
also the eigenvalue with largest real part of a deformation M{^) of the Markov matrix of the system [l^, (3]. The 
matrix M(pi) is similar to the hamiltonian of a XXZ spin chain with twisted boundary conditions [loj . and can thus 
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be diagonalizcd using the Bcthc Ansatz. The generating function of the cumulants of the current is given by 

^'^^'-(l-)4.o.f^) . (3) 

|t = l 



p ^ ' dt ^yx^-Qit/x) 



where the polynomial Q of degree n, along with a polynomial R of degree L, is a solution of the functional Bethe 
equation [17 1 



Q{t)R{t) = e^^il - t)^Q{xt) + (1 - xtfx'^Qitlx) . (4) 

This equation has several solutions, corresponding to different eigenstates of the matrix Af(7). The solution of 
equation ^ corresponding to the largest eigenvalue of M{j) verifies 

Q(t) = r + 0(7), (5) 

and 

e"'^Q(l) = x^Q{l/x) . (6) 
Using (O and the functional Bethe equation @ can be solved perturbatively near 7 = [isj . 

II. TREE STRUCTURES 

We now introduce a few combinatorial structures in terms of which the generating function of the cumulants of the 
current £^(7) will be expressed. We call "composite node" a (finite) set containing a strictly positive odd number of 
"elementary nodes". The size |c| of a composite node c is defined to be the number of elementary nodes it contains. In 
the following, two different sets of composite nodes will be considered: the set C of composite nodes without internal 
structure and the set C of composite nodes with an internal unrooted tree structure [3| (or acyclic graph structure) 
linking all the elementary nodes it contains. The elementary nodes and the composite nodes of size 1 will both be 
represented by a dot (•), the composite nodes elements of C\{#} by their size surrounded by a circle {e.g. (3), (5), . . . ), 
and the composite nodes elements of C\{#} by the tree structure on the elementary nodes they contain surrounded 
by a closed line {e.g. [% % •) ) . 

We can now build trees whose nodes will be composite nodes. We call Q the set of unrooted trees with nodes 
elements of C and Q the set of rooted trees with nodes elements of C and oriented edges labeled by ii, 12, ... One 
arbitrary elementary node of g G chosen as the root of 5, will be represented by a small circle (O). The edges 
between composite nodes are identified to edges between elementary nodes belonging to different composite nodes. 
These edges will be called "outer edges", by opposition to the "inner edges" linking elementary nodes contained in 
the same composite node. For g (z G or g G G, we define the size jgl of g as the sum of the sizes of the composite 
nodes of g. We call Gr (respectively Gr) the subset of G (resp. G) consisting of trees of size r. The first sets Gr and 
Gr are drawn in fig[l]and figlH 

For g £ G, call c{g) the set of the composite nodes of g. For a composite node c € c{g), we define Vc as the 



», (!>-•-•, •-(3>-», (5) 



FIG. 1: The five first Gr sets. 

number of composite nodes in g which arc neighbors of c. For g (z G, e{g) will be the set of the elementary nodes of 
g and o{g) the set of the outer edges of g. For an outer edge o S o{g), we call £{o) the label of o and m{o) the sum of 
the labels of g which label both the inner and outer edges of the subtree of g beginning at o, that is, the sum of all the 
labels that can be attained from o by moving on the edges of g away from the root of g. Finally, for an elementary 
node e of g € G, we will also use the notation £{e) for the sum of the labels of the edges pointing to e minus the sum 
of the labels of the edges leaving e. 
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FIG. 2; The four first Qr sets. The . . . represent trees obtained from the trees drawn by changing the root, the direction of the 
arrows on the edges and doing any permutation of the labels of the edges. 



For g G G, we define a tree g* G Q as follows: we attach all the elementary nodes of g to the root such that all the 
edges point away from the root and keep their labels. Then, we set the composite nodes such that all the composite 
nodes are of size 1 except the one containing the root, and such that the labels of the outer edges of g* arc the same 
as the labels of the outer edges of g. For example, wc have 



(7) 





For these trees, the values £{e) for the elements e of e{g) and e{g*) are respectively 



g ■ is, h - h, -h - «2, «2 - «4, ii + k - «6, -«5 - «7, ir, H + is, -is 
g* ■ ii, «2, «3, «4, ib, ie, i7, is, ~ii - «2 - «3 - «4 - «5 - ^6 - «7 - «8 , 

while the values m(o) for the elements o of o(g) and o[g*) are given by 

g : ig, Z7, ig, i4 + i5 + + ir + is • 
g* ■ iz, ii, ij, is ■ 



(8) 
(9) 



(10) 
(11) 



We will now define some functions acting on Q and Q. For g E G and two arbitrary functions if and r/, we define 



,eGe(g) / \e£e{g) 



(12) 



It can be shown that Uip^r/ig) docs not depend on the position of the root or on the composite nodes of g, but only 
on the tree structure of the elementary nodes of g. For an arbitrary even fmiction ^, wc also define 



vd9)= n ^Mo)) 

oeo(g) 



(13) 



Unlike U^^r){g), ^(ff) does depend on the position of the root and of the composite nodes in g, but not on the direction 
of the edges and on the internal tree structure of the composite nodes. We choose an arbitrary application from Q 
to Q preserving the tree structure on the composite nodes, which roots the trees, adds an internal tree structure on 
the composite nodes, labels the edges, and adds arrows to them. Then, we define for g G Gr 



(14) 



For r = 1, we take W^'^{9) = ip{0)ri{0). It can be proved that Wi^'^ does not depend on the choice of the function 6. 
Performing some changes of variables on the summation indices ii, . . . , v-i in equation p4p . an equivalent expression 
for W^'^{g) can be written: 



K''(9) = E • ■ • E u^AHgrmeig)) 



(15) 
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For a tree g G wc define a symmetry factor St{g) associated to g by 

wliere Pt(i?) is equal to the mimber of permutations of the composite nodes of g which leave it invariant. Table U gives 
a few examples of symmetry factors. These symmetry factors can be used to define a suitable generating function 
for the trees of Q. From a direct calculation up to r = 16, we observed that the following generating function can be 
expressed in a simple closed form: 

Zr{z) = E ^TT = ^ n [''(^ + 1) - ' (17) 
where \c[g)\ is the number of composite nodes in g. 

geg m •-• • • • (3) (5) ^)»-^ ^Ml) 

St{g) 12 2-18 200/3 8 -12 

Pt{g) 12 2 11 8 2 

TABLE I: Examples of symmetry factors of trees g £ Q 



III. PARAMETRIC EXPRESSION FOR THE CURRENT FLUCTUATIONS 

We now state our conjecture about the current fluctuations for the ASEP in terms of the trees and forests that we 
introduced before. We define the function ipi by 



the function jy by 



and the (even) function by 



, , {n + z) / L - n - z 
^'(")^ (L^T) I 



■q{z) = 



L 

n + z 



(18) 



(19) 



1 if z = 



(20) 



The value ^^^(O) was set to 1, but we emphasize that an arbitrary value could have been taken by modifying the 
following accordingly. With these definitions, the polynomial Q solution of (jlHS]) is then given by 



In 



Q{t) 



x-^Q(t/x) 



2^2^ 2k-i 2^ St(g) ' 



(21) 



with 



B = {-iy 



Q(0). 



(22) 



Using the method developed in [T^ to solve the functional Bethe equation Q perturbatively up to order 5 in 7 
and 1 — i, we checked equations (|2ip and (|22p at this order for all the systems up to size 12. Making the analytic 
continuation for complex L in equation (j2ip . the apparent divergences from ipi when L takes an integer value such 
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that L < I vanish since ipi{z) is always multiplied by rj^z). Taking t = 1 in (|2ip . the linear term of (po{z) = {n + z)/L 
does not contribute to W^'^^, and we obtain from ^ 

k=l ^ ^ geQk 

Taking the derivative at t = 1 of equation (|2ip . we obtain from ^ 

Eij)^J^ _ 2il-x) ^fB^^ W^^^^jg) 

where the mean value of the current J is given by J/p = (1 — x)n(L — n)/(i — 1). We used again the fact that the 
linear term of (fi cancels in W^f, while the constant term gives (I23p . 

Equations ([221) and ([M]) give a parametric expression for i?(7) similar to the one obtained for TASEP in [l^. In 
the TASEP limit a; = 0, ^^{z) = 1 for all z. Thus, from neither Wl^'^ig) nor W^i\g) depend on g G Gk anymore. 
They are equal to 



i^r(5)i..o = M;:) /c) (25) 



kn 



1 k{k-l)n{L-n) (kL\ (L 
W^. (5)|..o= ^^^^ U-J 



fc 



From the generating function of the trees Zk{z) at z ~ 1, we recover the known parametric expression [Tsj for 
i?(7) in the TASEP limit. We also recover the large L limit with non vanishing asymmetry considered in [l^ . 



IV. EXPLICIT EXPRESSION FOR THE CUMULANTS OF THE CURRENT 

The parametric expression (|23ll24p for the generating function i?(7) docs not give directly access to the cumulants 
of the current. To obtain the fc-th cumulant, one has to eliminate the parameter B between equations and (|24p . 
This can in fact be done systematically at all orders in 7. For this purpose, we have to introduce other combinatorial 
objects: forests, that is sets of trees. 

We call Ti. the set of forests with trees elements of Q of size strictly larger than 1. We call TL the set of forests with 
trees elements of Q of size strictly larger than 1 and with edges relabeled by zi, 12 ^ ... such that all the labels are 
different. The size \h\ of a forest h is defined to be the sum of the sizes of the trees it contains. The number of trees 
in a forest h will be denoted by h. We call Tir (respectively Tir) the subset of TL (resp. Ti.) with forests h such that 
\h\ — h = r. For h G TLr, the number of edges in h is equal to r. 

The various functions defined on Q and Q (e.g. the sets of composite nodes c, elementary nodes e, outer edges o, 
the application 9, . . . ) extend naturally to 7i and H. The operator * is defined for a forest h G TC as h* = {g*,g £ h}. 
The functions U^,rj, and VF^'^ are simply extended to forests h by replacing the e{g) and o{g) by the corresponding 
e{h) and o{h). Both expressions ([T^ and ([TS]) for W^'^ are still equivalent in the case of forests. We will also need a 
symmetry factor for forests. For h ^ H, we define 

where Pf{h) is the number of permutations of the identical trees in the forest h. 

The equation (|23p for 7 in terms of B can be inverted by using the Lagrange inversion formula (see 6.17. flst). The 
sum over trees in equation (|24p becomes then a sum over forests, and we finally obtain 

£;(7)^J7 _ 2(l~x) ^ (-L7r V ^'^^^ (07) 
p L(L-l)^ 2^r\ ^ Sf(h) ■ ^ ' 

^ ^ ' r=2 heHr-i ' ' 

This equation for £^(7) was checked by solving the functional equation ^ perturbatively (by the method described 
in [1^) up to order 7 in 7 for all the systems with 2 < L < 12 and n < L/2. The cases n > L/2 then follow from the 
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particle- hole symmetry of the system, which is verified by equation ([27]) . Equation ([27]) immediately gives the known 
results im for the three first cumulants, and also expressions for the higher cumulants. Up to order 4 in 7, we have 



48 



W 



+ 3W[m • • •] - [(3>-»] 



+ 0(7') , (28) 

with W = W'li-'' . 

In the large system size limit with finite density p = n/L, the expression (j27p for E{'^) simplifies. In the scaling 
1 — x ~ 2$/-\/Lp(l — p), Stirling's formula gives the leading order in the size L of the system for all the cumulants 
of the current: 



£(7) - J7 _ 4$ ^ (-7 - P)^') ' 



r2 Orr\ 2^ Rrih\ ' ^ ' 



p i2 2^r\ ^ Sf(h) 

Here, is defined in the same way as W except for the fact that the discrete sums over the indices ij are replaced 
by integrals between —00 and +00. 



V. CONCLUSION 



With our conjecture (j23ll24p for the generating function of the cumulants of the current, we recover the exact result 
by Derrida and Lcbowitz for the totally asymmetric model [l^ , as well as the expression obtained by Lee and Kim for 
the partially asymmetric model with non vanishing asymmetry [l6j . The expression (|27|) for E{j) then gives explicit 
expressions for the cumulants of the current, and we recover the known results for the three first cumulants [l5|. It 
would be interesting to calculate from (|27|) the cumulants of the current for the weakly asymmetric model, where the 
asymmetry 1 — x scales as 1/L. It should be possible to recover the expression obtained in [l^ in terms of Bernoulli 
numbers. So far, we only checked this for the four first cumulants. 

A distinctive feature between the known result for TASEP [l^ and the expression with partial asymmetry ([?7|) 
is the appearance of tree structures. This emphasizes the importance of combinatorics in theoretical physics, in 
particular in relation with integrable models. It would be interesting to know whether tree structures also appear in 
other situations, in particular for the open ASEP, the multispecies ASEP and the ASEP on the infinite line. 
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